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ABSTRACT 

We analyze the photometric redshift catalog of the Sloan Digital Sky Survey Data Release 5 (SDSS 
DR5) to estimate the Fisher information in the galaxy angular power spectrum with the help of 
the Rimes-Hamilton technique. It is found that the amount of Fisher information contained in the 
galaxy angular power spectrum is saturated at lensing multipole scale 300 < I < 2000 in the redshift 
range 0.1 <photo-z< 0.5. At I = 2000, the observed information is two orders of magnitude lower 
than the case of Gaussian fluctuations. This supports observationally that the translinear regime 
of the density power spectrum contains little independent information about the initial cosmological 
conditions, which is consistent with the numerical trend shown by Rimes-Hamilton. Our results also 
suggest that the Gaussian-noise description may not be valid in weak lensing measurements. 
Subject headings: cosmology: theory — methods:statistical — large-scale structure of universe 



1. INTRODUCTION 

Cosmic shear, which refers to the weak gravitational 
lensing by large-scale structure in the universe, is an 
invaluable tool in cosmology since it directly probes 
the gravitational clustering of dark matter in the uni- 
verse without concerns about light-to-matter bias. Very 
recently, the three dimensional map of dark matter 
distribution has been indeed reconstructed on cosmo- 
logical scales from the measurements of cosmic shear 
(jMassev et al.l 120071 ) . which marks an advent of lensing 
cosmology. Among many lensing observables, the angu- 
lar power spectrum is regarded particularly important, 
since it is in principle capable of constraining the key cos- 
mological parameters with high precisi on, including the 
eni gmatic nature of dark energy (see iHoekstra fc JainI 
20081 for a recent review). 

The success of the angular power spectrum as a cos- 
mological probe, however, is subject to one critical issue: 
how much information does it preserve about the initial 
conditions of the universe? A natural expectation is that 
some but not all information might have been destroyed 
in the subsequent nonlinear evolution. The amount of 
information contained in the power spectrum about the 
initial cosmological conditions is directly related to the 
precision of the cosmological parameters constrained by 
using the matter power spectrum. In other words, the 
non-Gaussian errors caused by the loss of information 
in the power spectrum would propagate into large un- 
certainties in the determination of the cosmological pa- 
rameters. The error propagation and precision on cos- 
mological parameters can be quantified in terms of the 
Fisher informa tion co ntained in the matter power spec- 
tr um (Tegmar k et al.l l l99 7). 

iRrines fc HamiltonI (f200a hereafter, RH05) have for 
the first time estimated the Fisher information in the 
matter power spectrum. Basically, RH05 have calculated 
the information content about the amplitude of the lin- 



ear power spectrum averaged over an ensemble of many 
N-body realizations and found very little independent 
information at the translinear scale. Somewhat surpris- 
ingly, however, RH05 also found that there is a sharp rise 
in t he amount of information in the nonlinear sc ale (see 
also iRimes fc HamiltonI 120061: iRimes et al.ll20Q6D . This 
phenomenon has been also noted by the analytic work of 
iNevrinck eTell (|2006D . 

Yet, according to the recent results derived by 
INevrinck fc Szapudil ()2007l ) based on the halo model, the 
Fisher information in the density power spectrum about 
all key cosmological parameters including the initial am- 
plitude is highly degenerate both in the translinear and 
the nonlinear regime. Their work has indicated that it 
might not be possible to extract the initial cosmological 
conditions from the nonlinear dark matter power spec- 
trum to a high statistical accuracy. 

Thus, the previous numerical and analytic results fore- 
cast that the accuracy in the determination of the cos- 
mological parameters in lensing cosmology may be lower 
due to the non-Gaussian errors. As the next generation 
of large surveys optimized for weak lensing will soon be 
on the pipeline, it is imperative to test observationally 
the information content in the lensing power spectrum. 
In this Letter we attempt to do this test by applying the 
RH05 technique to the photometric galaxy catalogs from 
the Sloan Digital Sky Survey Da ta Release 5 (SDSS DR5, 
lAdelman-McCarthv et aDl2007t l at typical lens redshifts 
z ~ 0.3. 

2. DATA 

We use the Photoz2 and the PhotoPrimary catalog 
which are both publicly available at the web site of the 
SDSSDR5 (http://www.sdss.org/dr5). The Photoz2- 
catalog contains information on galaxy photometric red- 
shift (photo-z, Zp), while the PhotoPrimary catalog has 
imaging data of galaxy right ascension (a), declination 
(Si), and dereddened model magnitude (Mr). For our 
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TABLE 1 



photo-z 


Si 


a 


# 


of galaxies 


# 


of cells 


0.01 


a 0.1SZp<0.2 

i o.aszp<o.3 

o 0.3gZp<0.4 

: » 0.4SZp<0.5 

with shot 


0.1,0.2) 
0.2,0.3) 
0.3,0.4) 
0.4,0.5) 


0°,60°) 
0°,60°) 
0°,60°) 
0°,60°) 


120°, 240°) 
120°, 240°) 
120°, 240°) 
120°, 240°) 




2794590 
4865953 
9768987 
9500664 




55 
55 
55 
55 



analysis, we select those galaxies in the catalog whose 
dereddened model magnitude, photo-z and angular posi- 
tions satisfy the constraints of Mr < 22.0; 0.0 < Zp < 1; 
0° <Si <60°; 120° <a<240°. 

Selected are a total of 42587179 galaxies which we di- 
vide into smaller samples according to their angular po- 
sitions and photo-z 's. We first decompose the selected 
portion of the sky into 55 cells, each having approxi- 
mately same linear size of 10°. To make equal-size cells 
on non-flat sky, the bin size of Sj is set at the constant 
value of Si = 10°, while the bin size of a is adjusted by a 
factor of 1/ cos Sim to compensate for the decrease in the 
width of Aa with Si. Here, we choose Sim as the median 
value of Si at each cell. We end up with having a total 
of 55 cells on the plane of the sky in the selected angu- 
lar range. Table [1] lists the ranges of a and Si, a total 
number of galaxies {Ng) and cells (Nc) in four photo-z 
slices. 

We treat approximately each cell as a square region 
and decompose it further into 128^ pixels of equal size. 
Counting the number of galaxies at each pixel, we find 
those pixels where no galaxy is located and mask them 
and their nearest neighbor pixels as gap regions, assign- 
ing zero density contrast. After masking the gap re- 
gions, we bin the selected photo-z range with bin size 
of Azp = 0.1. A set of 55 cells in each redshift slice rep- 
resents an ensemble of realizations for the calculation of 
the galaxy angular power spectrum. The size of the field 
is well in linear regime, so the small-size cells should be 
reasonably independent. Note that we don't include the 
effect of beat -coupling discussed in Rimes & Hamilton 
(|2006f ) and in lRimes et aH ()2006f ). 

For each cell in each photo-z slice, we construct the 
two dimensional density field on 128^ pixels by deter- 
mining the dimensionless density contrast S as residual 
number density of the galaxies, S = (n — n)/h, where n is 
the number density of the galaxies belonging to a given 
pixel, and n is the mean number density of galaxies av- 
eraged over a given cell in a given photo-z slice. For the 
calculation of the mean number density n, we exclude 
the gap regions, counting only unmasked pixels. 

Then, we perform the Fourier-transformation of the 
two dimensional density field with the help of the FFT 
routi ne provided in the numerical recipes (jPress et al.l 
1199^ and calculate the galaxy angular power spectrum, 
Cg{1), as a function of a multipole I in the Fourier space: 



Cg{1) 



-^total /ir/,\|2\ 



iV, 



unmask 



-mm, 



(1) 



where 1 is a two dimensional multipole vector represent- 
ing the Fourier transform of the two dimensional angu- 
lar position vector, iVtotai is the total number of pixels 
in a given ceU (i.e., iVtotai = 128^), and A^unmask repre- 
sents the number of pixels which do not belong to the 
masked gap regions in a given cell. We weight the an- 
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Fig. 1. — Mean angular power spectrum in a dimensionless form, 
C(l)l(l + l)/27r, averaged over 55 cells in four photo-z slices: 0.1 < 
2 < 0.2 (open squares); 0.2 < z < 0.3 (open triangles); 0.3 < 
z < 0.4 (open circles); 0.4 < z < 0.5 (open stars). The errors 
are calculated as one standard deviation in the measurement of 
the mean angular power spectrum. The dotted lines represent the 
cases before the shot-noise subtraction. 

gular power spectrum by a factor of (A^unmask/^^totai)""^ 
to compens ate for the zero powers assigned to the 
gap regions (jHamilton 1 120050 . The average fraction of 

A^unmask/A^total is fouud tO bc 0.92 ± 0.02. 

The galaxy angular power spectrum, Cg(Z), obtained 
by equation ([l} includes the shot-noise. To eliminate 
the shot-noise power from a given cell to which a total 
of Ng galaxies belong, we construct a random sample of 
equal number of galaxies by generating the values of (J/, 
a) randomly. If a randomly generated galaxy is found 
to fall in gap regions, we regenerate {5i,q) till it be- 
longs to a unmasked pixel. Using this random sample 
of Ng galaxies, we repeat the whole process: construct- 
ing the density field on 128^ pixels, Fourier-transforming 
the density field, calculating the shot-noise power spec- 
trum, Csn(0; which is also weighted by the same factor 
of A'totai /^unmask- Finally, we subtract the shot-noise 
power spectrum, Csn, from Cq{1), to obtain the shot- 
noise-free angular power spectrum, C{1) for all 55 cells 
in each photo-z slice. 

Figure [1] plots the mean angular power spectra, l{l + 
\)C{1)/2tt, averaged over 55 cells after the shot-noise 
subtraction as solid lines with errors for four different 
photo-z slices. The errors, crc, are calculated as one stan- 
dard deviation in the measurement of the mean angular 
power spectrum averaged over 55 cells belonging to each 
photo-z shce: ac = -sj {l^C^il)) /{Nc -I). For_the com- 
parison, the mean angular power spectrum, Cg{1), be- 
fore the shot-noise subtraction is also plotted as dotted 
lines. As can be seen, before the shot-noise subtraction, 
the slope of the mean angular power spectrum changes 
from '^ 1 to ~ 2 as the multipole I increases. After the 
shot-noise subtraction, its slope stays constant at ~ 1. 

The fluctuations of the angular power spectrum at the 
lowest multipole (/ < 100) represent the presence of the 
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sample variance caused by the finite size of the sky cells. 
Note that the amplitude of the angular power spectrum 
decreases with photo-z, indicating the growth of powers. 
Here, we focus mainly on the photo-z ranges of [0.1, 0.5): 
At higher photo-z {zp > 0.5), the shot-noise will domi- 
nate; At lower photo-z {zp < 0.1), due to the increase 
of the fractional distance errors, the three dimensional 
analys is will be necessary to measure the powe r spec- 
trum (jDodelson et al.ll2001l : [Tegmark et al.ll2002[ ). 

3. ANALYSIS 

The angular power spectrum estimated in §2 is sup- 
posed to contain information about the amplitude of the 
linear power spectrum: C{1) — C{l;\nA) where \nA 
is the log of the amplitude of the linear density power 
spectrum. According to RH05, the Fisher information 
contained in the angular power spectrum about the am- 
plitude of the li near power spectrum can be written as 
(|Tegmark et al.lll997| ) 



/ = - 



E 



dlnCjk) 
dliiA 



Cov7 



dlnA 



(2) 



Here (CoVj ) represents the inverse of the covariance 
matrix whose components are given as 



Cov; 



{ACik)AC{l,)) 
C{U)C[l,) ■ 



(3) 



where AC(I) = C{1) — C{1) is the scatter of the angular 
power spectrum from the me an value at a given multipol e 
/. As explained in detail bv iRimes fc HamiltonI (|2006D . 
for the case of linear power spectrum, C^iJ), from the 
Gaussian density fluctuations, the inverse of the covari- 
ance matrix is diagonal and equation ([3]) will be simpli- 
fied into I{li) — n{li)/2 since d\nC{l)/d\nA = 1 and 
Cov^ = n{li)/2 where n{li) is the number of Gaussian 
modes around the multipole k. 

RH05 asserted that the amount of information pre- 
served in the nonlinear power spectrum depends on the 
existence of an invertible one-to-one mapping between 
the linear and nonlinear regime. If such a mapping ex- 
ists, then the same amount of information that the lin- 
ear power spectrum contains about the initial ampli- 
tude can be extracted also from the nonlinear power 
spectrum. If not, then it would be difficult to deter- 
mine the initial amplitudes with high precision from the 
nonlinear power spectrum. However, note that what 
we have measured in §2 is the nonlinear galaxy power 
spectrum. Unlike the nonlinear matt er power spec- 
truni that is reasonably well u nderstood (jHamilton et al.l 
119911 : [Peacock fc Doddsl[T996[ ). the non-hnear clustering 
of galaxies is rather heuristic and usually constructed to 
fit the data. Although the non-linear galaxy power spec- 
trum C(/) is known to be neither the nonlinear matter 
power spectrum nor the linear matter power spectrum, 
we set C{1) to the linear matter power spectrum AC'^{1) 
for simplicity. 

Using the mean angular power spectra measured at 
four photo-z slices, we first calculate the normalized co- 
variance matrix defined as 



Gov,; 



{AC{h)AC{l,)) 
^{ACHk)){ACHl,)y 



(4) 



Figure [3] plots the magnitudes of (Gov^j) as grey scales 
for the cases of four photo-z slices: [0.1,0.2), [0.2,0.3), 
[0.3,0.4) and [0.4,0.5) in the panel (a), (b), (c) and (d), 
respectively. As can be seen, the magnitudes of the off- 
diagonal elements of (Govy) deviate from zero, which 
suggests the reduced amount of information than the case 
of Gaussian fluctuations. 

Finding the inverse of the covariance matrix. Gov" , 
we calculate the cumulative information /(< I) contained 
in the SDSS angular power spectrum: 



ii<h) = J2J2y^^''^^' 



Vj^ 



k^l. 



,ni 



(5) 



where y = {yi) is a column vector whose component is 
all unity and ni is the total number of multipole bins. 
Figure [H plots /(< I) for the cases of four photo-z slices. 
The errors are calculated from the 1000 times bootstrap 
resampling, and the dashed line represents the case of 
Gaussian- fluctuations. We offset slightly the horizontal 
positions of the four symbols at each multipole bin to 
show the errors clearly. 

If there were a one-to-one map from the nonlinear to 
the linear regimes, then J(< I) obtained from SDSS cat- 
alogs would follow the dashed line. As can be seen, at 
low multipoles {I < 300) it matches the Gaussian case 
pretty well. At higher multipole {I > 300), however, we 
detect a clear signal of saturation in /(< I). Note that at 
I = 2000 the observed amount of the cumulative infor- 
mation is two orders of magnitude lower than expected 
for the case of Gaussian-fluctuations in all four photo-z 
slices. The amplitude of the saturated region of /(< I) at 
high I tends to increase as photo-z increases in the range 
of 0.1 < Zp < 0.4, which is consistent with the numerical 
results of RH05. Although the amplitude becomes lower 
at the highest photo-z slice (0.4 < Zp < 0.5), the differ- 
ences are within the errors. Unlike RH05, however, we 
could not find a sharp rise of /(< I) at nonlinear scales 
even when we increase the number of pixels. It might be 
because in the nonlinear scale the light-to-matter bias be- 
comes important and many other complicated baryonic 
processes might have destroyed the information content 
further. 

4. DISCUSSION AND CONCLUSION 

It is worth mentioning here that what you have calcu- 
lated using the SDSS data is not rigorously the informa- 
tion about the amplitude of the power spectrum. But, 
it only approximates the information in the linear power 
spectrum to the extent that din P/ din A is near unity. 

Our result provides an observational evidence for the 
previous theoretical clues that at translinear scale there is 
very little independent information in the matter power 
spectr um, consiste nt with ligh t trac i ng matter till this 
scale (iRmTes fc Haniiito3[2005l[200l lRimes et aLll200l 
iNevrinck et al.l 120061 : iNevrinck fc Sza pudi 20d3). Since 
the multipole scales where the information saturation is 
detected correspond to the weak lensing regime, this re- 
sult has a direct impact on the weak lensing analyses. 
The loss of information in the angular power spectrum 
would lead to increasing non-Gaussianity in sampling er- 
rors. Unlike the usual assumption adopted in most weak 
lensing analyses that the non-Gaussianity contribution 
to the sampling errors for the angular power spectrum 
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Fig. 2. — Covariance matrices of the galaxy angular power spectrum, CU), at four photo-z slices: (a) 0.1 < Zp < 0.2, (b) 0.2 < 2; < 0.3, 
(c) 0.3 < Zp < 0.4 and (d) 0.4 < Zp < 0.5. As done in lRimes &: HamiltonI II2005I) . the grey scale is used to represent the magnitude of the 
covariances. 
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Fig. 3. — Cumulative information in the angular power spectrum 
for the four photo-z (zp) slices. The dashed line represents the 
cumulative information in the linear angular power spectrum and 
the errors arc calculated from 1000 times bootstrap-resampling. 



is marginal (| White k Hij[2000l : iCoorav fc Hul[200ll and 
references therein), our result suggests that it should be 
quite substantial and thus the Gaussian-noise description 
for the lensing power spectrum should not be valid. 
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